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ABSTRACT 

The relative distribution of abundances of refractory, intermediate, and volatile elements in stars with planets can be an important tool 
for investigating the internal migration of a giant planet. This migration can lead to the accretion of planetesimals and the selective 
enrichment of the star with these elements. We report on a spectroscopic determination of the atmospheric parameters and chemical 
abundances of the parent stars in transiting planets CoRoT-2b and CoRoT-4b. Adding data for CoRoT-3 and CoRoT-5 from the 
literature, we find a flat distribution of the relative abundances as a function of their condensation temperatures. For CoRoT-2, the 
relatively high lithium abundance and intensity of its Lii resonance line permit us to propose an age of 120 Myr, making this stars 
one of the youngest stars with planets to date. We introduce a new methodology to investigate a relation between the abundances 
of these stars and the internal migration of their planets. By simulating the internal migration of a planet in a disk formed only by 
planetesimals, we are able, for the first time, to separate the stellar fractions of refractory, intermediate, and volatile rich planetesimals 
accreting onto the central star. Intermediate and volatile element fractions enriching the star are similar and much larger than those of 
pure refractory ones. This result is opposite to what has been considered in the literature for the accreting self-enrichment processes of 
stars with planets. We also show that these results are highly dependent on the model adopted for the disk distribution regions in terms 
of refractory, intermediate, and also volatile elements and other parameters considered. We note however, that this self-enrichment 
mechanism is only efficient during the first 20 - 30 Myr or later in the lifetime of the disk when the surface convection layers of the 
central star for the first time attain its minimum size configuration. 

Key words. :planetary systems - Stars:abundances - Stars: fundamental parameters 

1 . Introduction with different properties from those of the Solar vicinity. Where 

does the bulk of the metal excess of the SWP mentioned above 



lGonzalezl ( ll997h and lSantos et"al](l2001h first identified a metal- originate? 

■ licity excess in stars with giant planets (SWP) relative to stars ^, ^ , , , ,. ,. 

» 1 » T-u » J u u c J • 1 The relative distribution of the stellar abundances of refrac- 

^ without planets. The tendency has been confirmed using large , • , , • , , ■ ,^ n • , 

i" I f,. XT »u 1 1 » u 1 u u J tory elements, with a high condensation temperature (rc), with 

samples of stars. Nevertheless, planets have also been observed , , •, , /, ^ , , ^, / 

^-s^ Ji .iiv • . J7=i — Z m respect to the volatile elements (low Tr), has been widely used 

around low metallicity main-sequence stars (ICochran et al.L .... , . . , „ , • , 

O" ItTwtK \TLr ^ u t u 1 4. 4. 1 t 1 in the literature as a tool to investigate the nature of metal enrich- 

. gOOTD. We note, however, that such a large statistical metal ex- c-nm / i ,'>f^f^-> 

__i . e A J f » ■ . »■ • » » -.u ment in SWP (see review papers of Gonzalez 2003, 200ob)). 

• cess found among dwarf stars is not present in giant stars with ^ ^ ^ y ^ jj 

^ planets (Pasquini et al 2007) enrichment mechanism should be triggered by the ac- 
• Using iron as a primary reference element, this metal excess, ^ion of a hot Jupiter The observed pile up of hot Jupiters is be- 
X has been the source of a number of studies trying to explain this Sieved to be the result of migration. These planets formed further 
b property. Two scenarios have been invoked: a self-enrichment ^^ay in the protoplanetary nebula and migrated afterwards to the 
■ mechanism of normal metal stars and a primordial scenario, in small orbital distances at which they are observed. During in- 
which SWP are the result of the formation of entire metal-rich ward migration, material (planetesimals) from the disk, depleted 
stars in equally metal-rich natal clouds. Which scenario is cor- ^ " ^e, could be accreted onto the star The presence of 
rect remains an open question to this day. ^ positive slope of [X/H] versus Tc, in which the abundances 
What about CoRoT exoplanet host stars? Although there are of ^ol^^ile elements are lower than that of the refractory abun- 
few cases with which to perform a statistical study, these stars dances, has been considered_a^ossible signature of a self- 
exhibit a different distributions: the majority have solar abun- enrichment mechanism (SmithetalJ,L2001). This is based on the 
dances, one has a mild metal excess, and one has a low metal "^^er P^^^n assumption that stellar accretion only favors near 
abundance. Even at this early stage, we are tempted to specu- ^^ar di sk regions ric h m refractory elements, 
late that CoRoT is exploring a different backyard of the Galaxy In [Winter et ap (|2007|), we explored by means of numerical 
simulations stellar accretion caused by an inward forced plane- 



Send offprint requests to: C. Chavero e-mail: carolinaOon . br tary migration in a disk formed only by planetesimals that cause 

* Based on observations made with the 2.2m telescope at the a metallic self-enrichment of the central star photosphere. At 
European Southern Observatory (La Silla, Chile). present, we attempt to identify this stellar accretion by deter- 
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mining the origin in the disk of these particles (planetesimals). 
In this way, we can separate the contributions of rocky ma- 
terial enriched differentially with refractory, intermediate, and 
volatile elements. In this work, we apply this methodology t o 



Table 1. Atomic parameters of the iron lines 



th ree different systems , namely : CoRoT-2 ( Alonso et aUl2008l) . 
3 jDeleuil etall l2008h . and 4 dAigrain et all l2008l) . The main 



purpose of this work is to analyze the mentioned slopes of 
[Fe/H] vs. Tc corresponding to CoRoT-2, CoRoT-3, CoRoT-4, 
and CoRoT-5 and their interpretations. For cases 2 and 4, we ob- 
tained new high resolution spectra. Cases 3 and 5 are taken from 
the literature. 

2. Observations 

The high-resolution spectra of CoRoT-2 and CoRoT-4 ana- 
lyzed in this work were obtained with the FEROS (Fiberfed 
Extended Range O ptical Spectrograph) echelle spectrograph 
dKaufer et all [1999) of the 2.2 m ESO telescope at La Silla 
(Chile), on October 23th, 2008. The FEROS spectral resolving 
power is R=48 000, corresponding to 2.2 pixel of 15 //m, and the 
wavelength coverage ranges from 3500 to 9200 A. The nominal 
signal-to-noise ratio (S/N) was evaluated by measuring the rms- 
flux fluctuation in selected continuum windows, and the typical 
values were S/N 55-60 after 2 x 3600 s of integration time. The 
spectra were reduced with the MIDAS pipeline reduction pack- 
age consisting of the following standard steps: CCD bias correc- 
tion, flat-fielding, spectrum extraction, wavelength calibration, 
correction of barycentric velocity, and spectrum rectification. 

3. Abundance analysis of CoRoT-2 and CoRoT-4 

3.1. Stellar parameters and iron abundances 

Determination of the basic, effective temperature (Tgff), surface 
gravity (log g), microturbulence (^,), metallicity ([Fe/H]), and 
the abundance analysis were derived by means of the standard 
approach of the local thermodynamic equili brium (LTE) u sing a 
revised v ersion (2002) o f the code MOOG (ISnedenL[T97 3) and 
a grid of iKuruca (Il993h ATLAS9 atmospheres, which include 
overshooting. 

The atmospheric parameters were obtained from the iron 
lines (Fe i and F n ) by iterating until the correlation coefficients 
between log e(Fe i) and lower excitation potential (xl), and be- 
tween log e(Fe i) and reduced equivalent width (log(W^//l) were 
zero, and the mean abundance given by F ei and Fe n lines were 
similar The iron lines taken from Lambert et alJ (11996) and 
ISantos et alJ (l2004l) were carefully chosen by verifying that each 
line was not too strong, and checked for possible blending. 

We adopted new \og gf values for the iron lines. These val- 
ues were computed from an inverted solar analysis using so- 
lar equivalent widths (EWs) measured from a solar spectrum 
taken with FEROS and a Kurucz grid model for the Sun (Kurucz 
1993) having {T^ff, log g, log e(Fe)) = (5777 K, 4.44 dex, 1 .00 
kms"', 7.47 dex). Table [1] contains the linelist used, the adopted 
parameters, and also includes the EWs measured in each star. 
Table |2] shows the atmospheric parameters derived by us and 
those found by other authors. 

We determined equivalent widths using the "splot" task of 
irafQ. Abundances were computed using the ABFIND driver 
in MOOG. 



Species 


A 


Xi 




EW(mA) 


EW (mA) 










CoRoT-2 


CoRoT-4 


Fei 


5044.22 


2.85 


-2.043 


88.0 


66.7 




5198.72 


2.22 


-2.176 




94.5 




5242.49 


3.63 


-1.130 


102.3 


84.1 




5322.05 


2.28 


-2.900 


77.2 


56.3 




5638.27 


4.22 


-0.809 


94.1 


77.6 




5731.77 


4.26 


-1.124 


73.0 


53.2 




5806.73 


4.61 


-0.893 


70.3 


51.9 




5852.22 


4.55 


-1.185 


52.2 


36.2 




5855.08 


4.61 


-1.530 


26.1 






5856.09 


4.29 


-1.568 


39.3 






6027.06 


4.08 


-1.190 


80.0 


61.6 




6056.01 


4.73 


-0.500 


92.2 


73.2 




6079.01 


4.65 


-1.010 




39.7 




6089.57 


5.02 


-0.884 


45.1 


36.5 




6151.62 


2.18 


-3.300 


70.6 


40.2 




6157.73 


4.07 


-1.240 


76.5 


61.3 




6159.38 


4.61 


-1.860 








6165.36 


4.14 


-1.500 


53.3 


42.0 




6180.21 


2.73 


-2.635 


71.7 


54.0 




6188.00 


3.94 


-1.635 


58.7 


45.6 




6200.32 


2.61 


-2.400 


98.0 


65.9 




6213.44 


2.22 


-2.592 


102.0 


76.2 




6226.74 


3.88 


-2.070 


33.6 






6229.24 


2.84 


-2.880 


45.5 






6265.14 


2.18 


-2.574 


111.0 


77.0 




6270.23 


2.86 


-2.588 


68.5 


49.0 




6380.75 


4.19 


-1.320 


66.9 


51.9 




6392.54 


2.28 


-3.933 


25.1 






6498.94 


0.96 


-4.632 


63.5 


38.0 




6627.55 


4.55 


-1.480 


32.1 






6703.57 


2.76 


-3.027 


43.1 


31.3 




6726.67 


4.61 


-1.050 


51.4 






6750.16 


2.42 


-2.621 


93.2 


75.4 




6752.71 


4.64 


-1.230 


47.0 


36.7 




6786.86 


4.19 


-1.900 


30.0 


21.0 


Fell 


5234.63 


3.22 


-2.238 


96.4 


100.8 




5991.38 


3.15 


-3.529 




47.9 




6084.11 


3.20 


-3.767 




33.3 




6149.25 


3.89 


-2.726 


44.9 


53.7 




6247.56 


3.89 


-2.342 


64.5 


72.2 




6369.46 


2.89 


-4.134 


23.7 






6416.93 


3.89 


-2.640 


43.1 


57.0 




6432.69 


2.89 


-3.564 


44.9 


52.9 




7711.73 


3.90 


-2.538 


49.8 





' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 



for Research in Astronomy, Inc.. 
National Science Foundation. 



under cooperative agreement with the 



Notes. Col 1. element; Col. 2: wavelength (in A ); Col. 3: excitation en- 
ergy of the lower energy level in the transition (in eV); Col. 4: oscillator 
strengths based on an inverse solar analysis.; Col. 5: EW for CoRoT-2; 
Col. 6: EW for CoRoT-4. 



3.2. Stellar abundances 

We present here results for 15 elements (O, Li, Na, Mg, Al, Si, S, 
Ca, Sc, Ti, Cr, Mn, Ni, Zn, and Ba), plus Fe. The abundances of 
Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Ni, Zn, and Ba were derived 
from the analysis of EWs for several unblended lines, measured 
by Gaussian fitting. To derive semi-empirical atomic \og gf val- 
ues for the lines of these elements, we used the EWs measured 
in the solar spectrum and performed an inverted solar analysis, 
as that carried out fort he iron element. The solar abunda nces of 
each element were taken from lAnders & Grevesse' (1989). Table 
[3] summarizes the linelist, the atomic parameters adopted, and 
the EWs measured in each star. For Li, O, and S, we used a 
spectral synthesis procedure. 
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Table 2. Stellar parameters of the CoRoT exoplanets parent stars 



Star Name 




logg 


^,[knVs] 


[Fe/H] 


mm 


V sin i [km/s] 


References 


CoRoT-2 


5696 + 70 


4.42 + 0.12 


1.71 ± 0.10 


0.03 ± 0.06 


0.03 ± 0.06 


8.5 ± 1 


This work 


CoRoT-4 


6115 + 70 


4.30 + 0.12 


1.37 ± 0.10 


0.12 ±0.05 


0.10 ±0.05 


5.5 ± 1 


This work 


CoRoT-2 


5625 + 120 


4.30 ± 0.20 






0.0 ± 0.1 


11.85 ±0.45 


Bouchv et al. (200S) 




5608 + 37 


4.71 +0.20 


1.49 ±0.06 


0.07 ± 0.04 






Ammler-von Eiff et al. (2009) 


CoRoT-4 


6190 ± 60 


4.41 ± 0.05 


0.94 ± 0.05 




0.05 ± 0.07 


6.4 ± 1 


Moutou et al. (2008) 


CoRoT-3 


6740 ± 140 


4.22 ± 0.07 




0.03 ± 0.06 


-0.02 ± 0.06 


17 ±1 


Deleuil et al. (2008) 


CoRoT-5 


6100 + 65 


4.189 ±0.03 


0.91 ±0.09 


-0.25 ± 0.06 


-0.25 ± 0.06 


1 ± 1 


Rauer et al. (2009) 



Notes. Column 1: star name; Col. 2: elfective temperature ; Col. 3: surface gravity ; Col. 4: microturbulence; Col. 5: metallicity; Col. 6: mean 
abundance (using element with more than three lines); Col. 7: projected rotational velocity Col. 8; references. 



It is important to study the highly volatile CNO elements be- 
cause due to their low Tq values, they are indicators of the cool 
disk zone below 371 K (Lodders, 2003). We took the conden- 
sation temperatures (Tq) from Lodders (2003i). covering a range 
between 70 and 1800 K. An attempt to measure abundances for 
an important volatile element, such as C, was unsuccessful. A re- 
liable calculation of the C abundance requires a higher S/N than 
used in this work. For instance, the C i lines in the region of 7 115 
A are hardly detectable and are severely blended. Similar prob- 
lems are found for other C i lines, but we instead measured other 
volatile elements as S and O for both stars. The abundance ratios 
[X/H] (and [Fe/H]) for each element, averaged over all useful 
lines, are presented in Table |4] along with the number of lines 
used in each case. 

3.2.1. Abundance uncertainties 

The internal eiTors in the adopted effective temperature and mi- 
croturbulence can be determined from the uncertainty in the 
slope of both the Fe i abundance versus excitation potential, and 
the Fei abundance versus reduced equivalent width relations. 
The standard deviation in logg was infeiTed by changing this 
parameter around the adopted solution until the and Fe ii mean 
abundances differ by exactly one standard deviation of the mean 
value of the Fe i abundance. These quantities are given in Table 

El 

Errors and uncertainties in the atmospheric parameters and 
determination for the continuum can affect abundance measure- 
ments in various ways. We tested the dependence of our results 
on atmospheric parameters obtained for CoRoT-2 and CoRoT4. 

Tables |5] and |6] show the variation in the abundance caused 
by A Teff, A log g, A and A W. The sixth column indicates 
the combined rms uncertainty of the second to fifth columns. 
The last column indicates the observed abundance dispersion for 
elements whose abundances were derived using three or more 
lines. 

We can observe that neutral elements are rather sensitive to 
temperature variations while single ionized elements are sen- 
sitive to the variations in logg. For the elements whose abun- 
dance is based on stronger lines, such as barium, the error intro- 
duced by the microturbulence is significant. For the elements an- 
alyzed by means of spectrum synthesis the same technique was 
used, varying T^ff, log g, and ^, then independently computing 
the abundance changes introduced by the variation in the above 
atmospheric parameters. 

The abundance uncertainties caused by the errors in the 
equivalent width measurements were computed fr om an expres- 
sion provided bv iCavrel de Strobel & Spita (Il988l) . The errors in 



the equivalent widths are determined, essentially, by the S/N and 
the spectral resolution. In our case, having ^ 48 000 and a typ- 
ical S IN of 55-60, the expected uncertainties in the equivalent 
widths are about 3 mA. These error estimates were applied to 
the measured W^s and the corresponding changes in the element 
abundances are listed in Col. 5 of Tables |5] and |6l 

The projected rotational velocities for both stars were deter- 
mined by applying the synthetic spectra method to several sin- 
gle and clean lines. The results presented in Table [2] are the fol- 
lowing: for CoRoT-2, we obtain vsin/ - 8.5 ± 1 kms ' and 
for CoRoT-4, we obtain vsin/ - 5.5 + 1 kms"'. These val- 
ues can be compared with t hose obtained for the same stars by 
other authors. For CoRoT-4. lMoutou et all (|2008[) found vsin / = 
6.4 + 1.0 kms w hich is simil ar to our value. The case of 
CoRoT-2 is different. iBouchv eT al. (2008) determined three in- 
dependent values for this parameter from different procedures: 
l)vsin/ = 9. 5 + 1.0 kms ' from SOPHIE cross-coiTelation func- 
tions (CCFs); 2) vsin J = 10.7 + 0.5 kms' from HARPS CCFs; 
and 3) vsin ; = 11.85 + 0.5 km s ' by modeling the radial veloc- 
ity anomaly occurring during the transit caused by the Rossiter- 
McLaughlin (RM) effect. The value obtained using SOPHIE are 
in closer agreement with our value since the difference is within 
the error bars. 

Column4 of Table|2]shows the values of the microturbulence 
calculated for each star. We found a small discrepancy between 
the f / that we f ound forCoRoT-4 (1.37 kms ') and that found by 
iMoutou et al] (12008) (0.94 kms '), although both works show 
an approximate solar value. This difference may be caused by 
the number of iron lines and the method used for determining 
this paramete r. ,Moutou et a l. (2008) used the automatic VWA 
method (Bru ntt et all 120081) . where they adjusted the microtur- 
bulence using an iterative procedure requiring that there is no 
correlation between the abundance and strength of a number of 
weak sensitive Fe I lines with equ ivalent widths < 80 mA. 

However, .Moutou et al. I (120081) are focused on presenting the 
discovery of the CoRoT-4b planet, and do not present a detailed 
analysis of the abundances values of the elements. We note, nev- 
ertheless, that our stellar parameters values, even if t hey have 
larger errors, are not incompatible with those of iMoutou et al.l 
(l2008h in terms of Teff and log g. 

3.3. Lithium 

Of the two observed stars, only CoRoT-2 exhibits a relatively 
strong lithium resonance line with a measured equivalent width 
of 124 mA. The secondary Lii line at 6104 A appears to be 
present in the spectrum as a very small feature in the wing of 
a Fe I line, confirming indirectly the high lithium abundance of 
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Table 3. Atomic parameters of the spectral lines used for each 
element. 



Species 


A 


Xi 




EW (mA) 
CoRoT-2 


EW (mA) 
CoRoT-4 


Nai 


5688.22 


2.10 


-0.678 


133.7 


110.3 




6154.23 


2.10 


-1.602 


43.0 


29.6 




6160.75 


2.10 


-1.316 


64.0 


43.0 


Mgi 


5711.09 


4.48 


-1.706 


120.0 


96.0 




6318.72 


5.21 


-1.996 




36.0 


Ah 


7835.31 


4.02 


-0.738 


43.0 


30.0 




7836.13 


4.02 


-0.569 


55.0 


50.0 


Si I 


5665.56 


4.92 


-1.980 


50.1 


33.0 




5948.55 


5.08 


-1.197 




81.8 




6142.49 


5.62 


-1.518 


35.3 


28.0 




6145.02 


5.62 


-1.420 


41.0 






6721.84 


5.86 


-1.152 


52.8 


48.0 


Cai 


5581.97 


2.52 


-0.728 


122.8 


92.0 




5590.12 


2.52 


-0.815 




82.0 




5867.56 


2.93 


-1.599 


33.0 






6166.44 


2.52 


-1.171 


84.0 


63.0 




6169.05 


2.52 


-0.812 


113.6 


90.0 




6169.56 


2.52 


-0.567 


145.7 


103.0 




6455.60 


2.52 


-1.420 


72.1 


45.0 




6471.66 


2.53 


-0.848 


118.0 


86.0 


Sen 


5526.82 


1.77 


0.111 


85.0 


83.0 




6245.65 


1.51 


-1.051 


37.0 


33.0 




6604.60 


1.36 


-1.146 




41.5 


Til 


5219.70 


0.02 


-2.268 


46.4 






5299.99 


1.05 


-1.425 


29.4 


16.0 




5866.46 


1.07 


-0.827 


70.0 


34.0 




6126.22 


1.07 


-1.398 


32.7 


12.0 




6258.11 


1.44 


-0.450 


59.9 




Cri 


5296.70 


0.98 


-1.411 


115.0 


84.0 




5300.75 


0.98 


-2.149 


70.0 


51.0 




5345.81 


1.00 


-1.036 




106.5 




5348.33 


1.00 


-1.301 




88.0 




5783.87 


3.32 


-0.220 


60.0 


40.0 




5787.93 


3.32 


-0.178 


55.0 


43.0 


Mm 


4265.92 


2.94 


-0.442 




52.6 




4470.13 


2.94 


-0.585 




43.0 




5399.47 


3.85 


-0.098 


44.2 






5432.54 


0.00 


-3.636 


60.0 


33.0 


Nil 


5084.11 


3.68 


-0.169 


106.0 


86.8 




6327.60 


1.68 


-3.100 


43.2 


38.1 




6643.64 


1.68 


-2.020 


111.0 


92.3 




6767.77 


1.83 


-2.170 


91.0 


77.5 




6772.32 


3.66 


-0.959 


54.8 


49.4 




7788.93 


1.95 


-1.951 


101.2 


88.3 


Zni 


4722.16 


4.03 


-0.390 


68.0 


70.0 




4810.53 


4.08 


-0.310 


87.0 


75.0 


Ban 


5853.69 


0.60 


-0.848 


90.5 


75.0 




6141.73 


0.70 


0.146 


160.0 


130.0 




6496.90 


0.60 


-0.222 


136.0 


109.0 



Notes. Column 1: element; Col. 2: wavelength (in A ); Col. 3: excita- 
tion energy of the lower energy level in the transition (in eV); Col. 4: 
oscillator strengths based on an inverse solar analysis.; Col. 5: EW for 
CoRoT-2; Col. 6: EW for CoRoT-4 



this star The lithium abundances were derived from synthetic 
spectra matching the observed Li i 6708 A resonance doublet. 

The wavelengths and oscillator strengths for the individ- 
ual hyperfine and isotope components of the lithium lines were 
taken from Smith et al. (1998) and Hobbs et al. (1999). The 
CN lines in the vicinity of the Lii doublet were included in 
the linehst. The solar "'Li/^Li isotopic ratio (''Li/^Li= 0.081) 



Table 4. Abundances of CoRoT-2 and CoRoT-4 



Species 


CoRoT-2 


#of 


CoRoT-4 


#of 




[X/H] 


lines 


[X/H] 


lines 


Oi 


-0.13 ± 0.10 


3 


-0.03 ± 0.10 


3 


Nai 


0.02 ± 0.03 


3 


0.05 ± 0.10 


3 


Mgi 


0.04 


1 


0.1 


2 


Ah 


-0.05 


2 


-0.02 


2 


Si I 


0.04 ± 0.05 


4 


0.01 + 0.08 


4 


Si 


0.00 ± 0.09 


1 


0.15 ± 0.08 


1 


Cai 


0.1 ± 0.05 


7 


0.1 ± 0.05 


7 


Sen 


-0.04 


2 


-0.03 ± 0.05 


3 


il I 


0.09 ± 0.09 


J 


O.OO + U. 11 


3 


Cri 


0.01 ± 0.08 


4 


0.1 ±0.04 


6 


Mm 


-0.06 


2 


0.02 ± 0.03 


3 


Fez 


0.03 ± 0.05 


32 


0.12 ±0.05 


27 


Fen 


0.02 ± 0.07 


7 


0.12 ±0.05 


7 


Nil 


-0.08 ± 0.06 


6 


0.17 ±0.04 


6 


Zni 


-0.1 


2 


0.01- 


2 


Ban 


0.19 ±0.04 


3 


0.14 ±0.06 


3 



Table 5. Abundance uncertainties of CoRoT-2. 



Species ATefj Alog g AW,, (S o"") CTobs 

+70 K +0.12 +0.10 3 mA 

Fei 
Fen 
Nai 
Mgi 
Ah 
Si I 
Cai 
Sen 
Til 
Cri 
Mm 
Nil 
Zni 

Ban 

Lii 006 OOO Mi ~ 006 005 

Oi -O03 O07 OOO - 0.08 0.10 

Si -0.05 0.05 0.00 - 0.07 0.09 



Notes. The second column indicates the variation in the abundance 
caused by the variation in Teff . The other columns refer, respectively, to 
the variations caused by to log g, f,, and W,i. The sixth column indicates 
the combined rms uncertainty of the second to fifth columns. The last 
column indicates the observed abundance dispersion for those elements 
whose abundances were derived using more than two lines. 

(lAnders & Grevess5 Il989h was adopted in the synthetic spec- 
trum calculations. 

For CoRoT-2, we obtained log e(Li) = 2.70 + 0.05. The ob- 
served and synthetic spectra of CoRoT-2 in the region of Li i 
line are shown in FiglT] For the secondary Li i line, we only 
obtained the upper limit to the Li abundance of log e(Li)< 2.8, 
which agrees with that obtained for the resonance line. The high 
Li abundance determined for CoRoT-2 of effective temperature 
Teft - 5696 K and both its chromospheric activity and relatively 
high rotation velocity vsin/ = 8.5 kms"' imply a young age of 
this star 

The star with the higher effective temperature reff= 61 15 K, 
CoRoT-4, exhibits a relatively weak Lii resonance line. 
Analyzing the relatively noisy spectrum of this star in Fig|2l we 
were unable to obtain a reliable value of the Li abundance. For 
CoRoT-4, we estimate only an upper limit to the lithium abun- 
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-OOl 
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0.04 


-0.02 


-OOl 


0.04 
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0.02 


-0.02 


-OOl 


OOl 


OOl 


OO 
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-0.02 
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0.0 


0.05 


-0.02 


0.08 


0.0 


-OOl 
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-OOl 


-0.03 


0.07 


OO 
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0.05 


OO 


-0.03 


OO 


0.02 


-0.03 


0.02 


0.03 


-0.05 
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0.05 


0.06 


0.08 


0.07 


0.03 


0.05 


0.03 


0.03 


0.06 




0.03 


0.04 




0.04 


0.04 


0.05 


0.04 


0.07 


0.05 


0.05 


0.07 




0.05 


0.09 


0.09 


0.05 


0.08 


0.08 


0.05 


0.09 




0.05 


0.08 


0.06 


0.05 


0.06 




0.05 


0.08 


0.04 
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Table 6. Abundance uncertainties of CoRoT-4. 



Species 


A Teff 


A l02 2 


Af 


A W,i 


) 


O" b 




_i_7n 1^ 
+ /U Jv 


_i_n 1 9 


_i_n 1 








Fei 


0.05 


-0.01 


-0.03 


0.05 


0.08 


0.05 


Fell 


-0.01 


0.04 


-0.04 


0.05 


0.08 


0.05 


Nai 


0.04 


-0.01 


0.0 


0.05 


0.06 


0.10 


Mgi 


0.04 


-0.01 


-0.01 


0.04 


0.06 




Ah 


0.03 


-0.01 


0.0 


0.05 


0.06 




Si I 


0.03 


0.0 


-0.01 


0.05 


0.06 


0.08 


Cai 


0.04 


-0.02 


-0.03 


0.04 


0.07 


0.05 


Sen 


0.0 


0.04 


-0.03 


0.05 


0.07 


0.05 


Til 


0.06 


-0.01 


-0.01 


0.08 


0.10 


0.11 


Cri 


0.06 


-0.01 


-0.04 


0.06 


0.09 


0.04 


Mm 


0.06 


0.0 


-0.01 


0.06 


0.09 


0.03 


Nil 


0.06 


-0.01 


-0.04 


0.05 


0.09 


0.04 


Zni 


0.03 


0.0 


-0.05 


0.06 


0.08 




Ban 


0.02 


0.02 


-0.08 


0.05 


0.10 


0.06 


Li I 


0.05 


0.00 


0.00 




0.05 




Oi 


0.12 


0.07 


0.08 




0.16 


0.10 


Si 


-0.03 


0.05 


0.00 




0.06 


0.08 



Notes. The second column indicates the variation of the abundance 
caused by the variation in Teff. The other columns refer, respectively, to 
the variations caused by logg, f,, and W,i- The sixth column indicates 
the combined rms uncertainty of the second to fifth columns. The last 
column indicates the observed abundance dispersion for those elements 
whose abundances were derived using more than two lines. 



dance on the order of loge(Li) < 2.2. This apparently high Li 
abundance of CoRoT-4 is, however, consistent with the Li abun- 
dance of the main-sequence field and cluster stars with effective 
temperatures of 5900 < Tf^f f < 6300 K (the so-called Li-plateau 
region according to .Boesgaard et al ] (120041 ) 

3.4. Oxygen 

To determine of oxygen abundance, we used the O i triplet lines 
at 7771, 7774, and 7 775 A . The at omic data were taken from 
the VALD database dPiskunov et al.L fl995; Kupkaet al., 19991). 
The value of the oscillation strength of the main component was 
modified to fit the solar abundance of loge(O) = 8.93 (Anders 
& Grevesse 1989). 

The oxygen triplet lines are known to be affected by NLTE 
effects, which lead to an overestimation of the oxygen abun- 
dances. Even if the forbidden [O i] line at 6300.304 A is assumed 
to be free from NLTE effects, this line is blended with a telluric 
[O i] emission line in our spectra and cannot be used for the oxy- 
gen abundance determinatio ns. To account for NLTE effects, we 
used the theoretical work of lTakedal (l2003i) . who calculated the 
values of the NLTE connections for each line of the O i infrared 
triplet for a grid of different atmospheric parameters and equiv- 
alent widths of the oxygen lines. 

For the cooler star of our sample, CoRoT-2, the NLTE cor- 
rection is about 0. 1 dex, whereas for the hotter star, CoRoT-4, 
NLTE effects are greater, about 0.2 dex. By applying the NTLE 
corrections, we derived the oxygen abundances log e(0) - 8.80 
and 8.90 for CoRoT-2 and CoRoT-4, respectively. 

For CoRoT-2, not only a strong resonance Li i line but also 
emission in the Ca ii lines is a strong indication of its youth and 
surface activity (see Fig. O. In the case of CoRoT-4, emission 
is absent (Fig. |3]l. This age-related activity is proposed to af- 
fect the oxygen abundances derived from the oxygen triplet lines 



m CoRoT exoplanet host stars 




6700 6702 6704 6706 6708 6710 6712 
WavelengT.h, A 

Fig. 1. CoRoT-2: Observed (dotted red line) and synthetic (solid 
blue lines) spectra in the region of the Li 1 6708 A line. Synthetic 
spectra were used to calculated the lithium abundances (from top 
to bottom): log e(Li) = 2.6, 2.7, and 2.8,. The best-fit function to 
the observations is achieved for lithium abundance log e(Li) = 
2.7 + 0.05. We also show three iron absorption lines. 
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Fig. 2. CoRoT-4: Observed (dotted red line) and synthetic (solid 
blue lines) spectra in the region of the Li 1 6708 A line. Synthetic 
spectra were used to calculated the lithium abundances (from 
top to bottom): loge(Li) - 2.1, 2.2, 2.3 and 2.4. It is difficult to 
determine a value for the abundance because of the noisy spec- 
trum, and we only estimated an upper limit of log£(Li) = 2.2. 
We also show three iron absorption lines. 



Chavero et al.: Abundance distributions in CoRoT exoplanet host stars 



Ca II (K) 

V 



mm 
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MB 
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Fig. 3. Calcium II K and H lines for stars CoRoT-2 and CoRoT- 
4. The central emission features in the Ca lines correspond to the 
younger and chromospheric active CoRoT-2 star. 



(IShen et all l2007h and may also affect the precision of the oxy- 
gen abundance determination for CoRoT-2 star. 

3.5. Sulfur 

Sulfur is a very important volatile element, being one of the 
commonly and observed elements, the only reference element 
i n an extended interval of Tc values between 704 K and 37 1 K 
dLoddersl l2003l) . To obtain the S abundances in both stars, we 
selected the line of S i at 6757 A. This line consists of three fine 
structure components and has to be studied using the syntheti c 
spectrum method. We used the \og gf data of lWiese et al.l(ll969h . 
which were provided by VALD. 

The value of the oscillation strength of the main compo- 
nent (-0.310) was modified to -0.340 to fit the solar abun- 
dance of loge(S) = 7.21 (Anders & Grevesse 1989). The sul- 
fur abundances that we derived are loge(S) - 7.21 + 0.09 and 
log e(S) = 7.36 + 0.08 for CoRoT-2 and CoRoT-4, respectively. 



4. Numerical simulationsianalysis and 
methodology 

4. 1 . The planetesimal disk 

A protoplanetary disk contains a mixture of gas and dust grains. 
These dust grains collide and build up larger and larger bodies, 
and small solid bodies known as planetesimals are presumably 
formed. Protoplanets grow within protoplanetary disks via pair- 
wise accretion of planetesimals. When the mass of a protoplanet 
is of the order of 10-20 terrestrial masses, it can gravitation- 
ally capture gas-rich disk, thus producing Jovian-type planets. 
Protoplanetary disk is expected to lose its gas component in the 
first 10 Myr, leaving a disk formed essentially of dust and plan- 
etesimals. 



For our purposes, we assumed that large planets observed 
close to their parent stars actually formed at larger distances 
but migrated inward due to energy loss and excess energy used 
to disperse the planetesimals (iMurrav et all Il998l) . A detailed 
descripti on of the planet mi gration used in this work is pre- 
sented in I Winter et al.l (l2007h . It consist of a forced inward mi- 
gration of a planet supposedly interacting with a planetesimal 
disk. Several possibilities were considered in terms of the planet 
eccer itricities, migration rates, and mass. A more general result 
from IWinter et al.l (12007 ) is that 2:1 mean motion resonance is 
the main mechanism for driving planetesimals to the surface 
of the star. This mechanism operates for slow migration rates 
(r - 10^ - 10^ yr) and low planet eccentricities. 

To construct a rocky disk model valid for ages equal to or 
older than 10 Myr, several considerations must be taken into ac- 
count. This disk containing a planetesimal distribution can be 
modeled using, for instan ce, knowledge obtained about a solar 
primeval nebula. Based on'LoddersI (|2003|) . in TableQwe present 
the distribution of the most important and more abundant rock- 
forming elements. We indicate the type of elements (from highly 
refractory up to highly volatile), their proper elements, their in- 
terval of condensation temperature, and finally, their semi-major 
axis distribution (in AU) with respect to a central solar type star 
whose radial thermal distribution is valid for an age older than 
10 Myr. 

For simplicity, we assume that in a certain epoch of the disk 
evolution, most probably in the range 20 - 30 Myr (see Sect. 5), 
the disk is formed by a sea of planetesimals defining three repre- 
sentative zones: a refractory rich planetesimals zone R between 
~ 0.03 AU and ~ 0.1 AU (between 1780 K and 1360 K), an in- 
termediate zone I for particles between ~ 0. 1 AU and 1 .56 AU 
(1360 K to 200 K), and a zone V of volatiles between 1.56 AU 
and ~ 4.5 AU ( < 200 K). 

In reality, the situation must be far more complex. Because 
nebular turbulence and dust migration probably occur during 
the earlier nebular stages, a radial mixture is expected. If zone 
R contains the pure refractory elements, zone / must contain a 
true mixture of elements. An example of this mixture in zone 
/, is the presence of Triolite (FeS) at 0.15 AU with a Tq 700 
K (^.odders, 2003). At this distance, an important condensation 
of pure Fe grains and a mild c ondensation of pure C may also 
occur ( Wehrstedt & GailLl200l . The zone / also contains the in- 
ternal end of the migration condensation front of H2O p articles 
at ~ 0.8 AU coming from the external zones of the disk ( IDavisl 
2005). In summary, we can consider that zone / contains plan- 
etesimals composed of a mixture of refractory type grains and 
volatile ices. The zone V is characterized by the abundant highly 
volatile CNO elements for which the highest Tq is only 180 K 
for oxygen. 



Table 7. Most important and more abundant rock-forming ele- 
ments (from Lodders (2003)) 



Type of element 


Element 


TcCK) 


a(AU) 


Refractory 
Common 

Moderately Volatile 

Volatile 

Highly Volatile 


Ca, Al, Ti 
Fe, Si, Mg 

P 

S 

H, C, N, + 
He, Ne, Ar 


1650 - 1360 
1360 - 1290 
1290 - 704 
704 - 371 

<371 


0.03 - 0.04 
0.04 - 0.046 
0.046-0.15 
0.15-0.48 
0.48 - < 5 
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4.2. The planet migration 

We numerically integrated the CoRoT-2, CoRoT-3, and CoRoT- 
4 systems within the circular restricted three-body problem of 
star-planet-planetesimals. The orbital and physical parameters 
for these stars are listed in Table [8] Each parameter is listed 
along with its formal uncertainty. We assumed that the planet 
may have been formed at a distance between 2 and 5 AU, where 
they are supposed to form under conditions where condensible 
materials exist. Thus, the planet is set initially with semi-major 
axes ap, =2, 3, 4, and 5 AU and eccentricity ep — 0.0. The planet 
is forced to migrate inward to its current position, ap^ , with con- 
stant speed in a time-scale r. We assumed that the planetesimal 
disk is composed of 1000 test massless particles moving around 
the star on initially randomly distributed circular orbits. The disk 
is interior to the orbital radius of the planet. The newly formed 
plan et is surrounded by an annular gap in the planetesimal disk 
(IWisdom. 1980; Lissauer, 1993). We assumed here an annular 
gap of 0.5 AU relative to the outer disc edge. Therefore, the outer 
radius is, in all cases, 0.5 AU smaller than the initial position of 
the planet. Thus the disk outer radius (ao/) is located at 1.5, 2.5, 
3.5, and 4.5 AU, according to the planet initial semi-major axis. 
The inner radius (ao,) of the disk is located at 0.03 AU. 

We are interested in studying the relationship between the 
metallicity of these CoRoT systems and the migration result- 
ing from interaction with planetesimals. Since a part of the mass 
of planetesimals will be accreted by the star during migration, 
we need to estimate for these CoRoT systems the mass of the 
disk contained inside the orbit of each planet. Considering the 
initial semi-major axes of the planets, we calculated the orders 
of the total mass Mq of the planetesimal disk required for the 
planet to migrate to its current position (Adams & Laughlin, 
|2003)- These data are presented in Table |9] We assumed that 
about 1% of the disk is in the form of solid bodies called plan- 
etesimals (iDavisL 12005 1). We note that the order of the total 
mass of the planetesimal disk increases with planetary mass and 
larger ap.. Using these values, we now estimate the percentage 
of the element types that fall on the star during the planet mi- 
gration. The numerical code SWIFT dUevison &Duncani[T994h 
was used. The integrator is based on the MVS method developed 
by Wisdom & Holman ( 1991). 

In general, the numerical simulations were performed for a 
timescale of 10^ years. The integration for each particle was in- 
terrupted whenever one of the following situations occurred: 

(i) collision between the planetesimal and the planet, i.e., 
when the planetesimal becomes closer than 1.2 planetary radii 
from the planet; 

(ii) collision between the planetesimal and the star, when the 
planetesimal becomes closer than 2.5 stellar radii from the star; 

(iii) ejection from the system, when the planetesimal reaches 
more than 50 AU from the star. 

4.3. Results 

The results of the numerical simulations are presented in Fig. 
m In that figure, we show the dependence of the percentage of 
planetesimals that fall on the star in terms of the value of a p.. The 
main significance of the results are the statistics in terms of per- 
centage of collisions with the star of planetesimals characterized 
by R, I, and V as a function ap.. The values of the percentages 
were obtained by measuring the number of planetesimals of the 
R, I, and V zones that collided with the star. Thus the percentage 
of collisions of planetesimals with the star gives an estimate of 
the mass of each type of material that the star will accrete due 



to the migrating planet. We note that while this percentage pro- 
vide reasonable estimates of the accreted mass, the assumption 
of there being a constant migration is not strictly correct since 
these collisions will tend to decrease as the planet migrates in- 
wards. 

We can also note from Fig.|4]that the percentage of planetes- 
imals in the volatile zone that fall on the star is increases with 
increasing value of ap.. This is due to the disk model considered. 
As mentioned in the previous section, the circumstellar disk is 
assumed to contain highly volatile particles between 1.56 and 
4.5 AU. They are distributed in that region for values of the ini- 
tial position of the planet beyond 2 AU. Thus, there are fewer 
particles in the volatile zone when the planet has an initial posi- 
tion 2 AU. In this case, all the particles are practically initially 
distributed in refractory and intermediate zones. Since the size of 
the refractory zone is much smaller (0.03 -0.1 AU) than the size 
of the intermediate zone (0.1 - 1.56 AU), there is a much larger 
quantity of particles in the intermediate zone. A more signifi- 
cant number of particles is then expected to fall on the star from 
the intermediate zone. Therefore, our results show that for these 
CoRoT systems, particles from the intermediate zone dominate 
the collision process onto the star In total about 79% of plan- 
etesimals from the intermediate zone were found to be falling 
onto the star In the CoRoT-2 and CoRoT-3 cases, all material (~ 
5%) of planetesimals from the refractory zone fell onto the star. 

For all the migration cases of the planet of the CoRoT-4 sys- 
tem, we note that planetesimals from the refractory zone did not 
fall onto the star. This occurred because the planet migration 
halts before the 2:1 resonance reaches the refractory zone. 

When planets migrated from the initial position at 3 AU, 
about 3% of the planetesimals from the refractory zone fell onto 
CoRoT-2 and CoRoT-3 stars. Planetesimals from the volatile and 
intermediate zones fell on the stars. The results show that about 
13% belongs to the volatile zone, while 59% belongs to the in- 
termediate zone. For the initial positions of the planet at 4 and 5 
AU, our results do not show significant differences between the 
values of percentages of planetesimals from the refractory, inter- 
mediate, and volatile zones that fall on CoRoT-2 and CoRoT-3. 
In these two cases, the numerical estimates are very close. There 
is doubtless no dependence of these percentages values with the 
initial semi-major axis of the planet. In general, giant planets are, 
therefore, be expected to form at distances between 4 and 5 AU. 
We focused on the situation of an orbital radius at 5 AU, where 
giant planets form and gain most of their mass. 

Motivated by the much smaller total mass of its planetesimal 
disk than the other systems (Table |9]l, the numerical simulations 
for the CoRoT-4 case were also performed for timescales of 10^ 
and 10^ years. For very rapid migrations, the planetesimals ex- 
perience close approaches to the planet and a significant nu mber 
of them are ejected from the system (Wi nter et al.L l2007h . The 
question that we wish to answer is: how much mass from the disk 
will the star accrete during these times? The results are shown in 
Fig- 13 In this case, we can note that 26% of planetesimals from 
the intermediate zone collided with star for t = 10"^ years, while 
for T - 10^ years the value is about 3.9%. 

4.4. Parameters dependence 

Until now, the results presented are concerning a given set of 
parameters model. In order to gain more insight into the numer- 
ical results from our disk simple model, it is useful to character- 
ize and, when possible, quantify its dependence with respect to 
some of the parameters used. To investigate the sensitivity of our 
results, in this section we present a study about the influence of 



7 



Chavero et al.: Abundance distributions in CoRoT exoplanet host stars 
Table 8. Orbital and physical parameters of the stellar systems 



Orbital parameter 


M,„„, (Mq) 


M,,;„„(,, (Mjupiter) 


Semi-major axis (AU) 


Eccentricity 


References 


Corot-2b 
Corot-3b 
Corot-4b 


0.97 ± 0.06 
1.37 ±0.09 

1 +0.03 


3.31 ±0.016 
21.66 ± 1.0 
0.72 ± 0.08 


0.0281 ± 0.0009 
0.057 ± 0.003 
0.09 ± 0.001 


0.0 (fixed) 
0.0 (fixed) 
O-O-co^"' 


Bouchv et al. (2008) 
Deleuil et al. (2008) 
Moutou et al. (2008) 



Table 9. Total mass of planetesimals for the planet migration up to its current position in the stellar system 



Disk mass (Mu) 


Corot-2b 


Corot-3b 


Corot-4b 


Up. (AU) 


2-5 


2-5 


2-5 


Mohk (Mq) 


0.013-0.016 


0.075 - 0.093 


0.002 - 0.003 




-I ' \ ' • ' ■ ' m ' m u-i , ^ , « : ■ : m ■ m , ^ , « , « , « , — 

012345 012345 012345 

% \ \ 



Fig. 4. These figures show the percentage of planetesimals that collided with the star of the CoRoT systems. Each curve gives the 
percentage of planetesimals of the zones refractory (square), intermediate (circle) and volatile (triangle). The values are given as a 
function of the initial semi-major axis of the planet, ap.. At the top of each figure, we indicate the corresponding system. In the first 
column, the results are for the case of zones R between 0.03 AU and 0.1 AU, I between 0.1 AU and 1.56 AU, and V between 1.56 
AU and 4.5 AU. The second column is for the case of zones R between 0.03 AU and 0.3 AU, I between 0.3 AU and 1.56 AU, and 
V between 1.56 AU and 4.5 AU. The third column is for the case of zones R between 0.03 AU and 0.1 AU, I between 0.1 AU and 
1.2 AU, and V between 1.2 AU and 4.5 AU. 



some of the adopted parameters on the amount of planetesimals 
that fall on the star. 



4.4.1. R, V and, I location zones 

A starting point for our discussion is provided by the distribution 
model of planetesimals in the disk, which is given by the location 
of the R, I, and V zones. This has important consequences. It im- 



8 



Chavero et al.: Abundance distributions in CoRoT exoplanet host stars 



CoRoT-E3:o4 



CoRo T-E:50 4 



CoRoT-Bso 4 



■ Refractary 

■ Imerni eiati 
H Volatile 



Type ti elements 

T„,^=100 000 yr 



Type ff elementB 

T„,j,=10 000 yr 



Type if elementB 

=1000 yr 



Fig. 5. Histogram of the percentage of planetesimals from the refractory (R), intermediate (I), and volatile (V) zones that fall on Corot-4 for the 
case of cip. = 5 AU. 



plies that the percentage of the type of elements that fall on the 
star seems to be uniquely controlled by the location zones con- 
sidered. As mentioned in Section 4.1, in reality a radial mixture 
of bodies of zones R and V must be contained in the whole disk. 
The intermediate zone is characterized by this mixture. Thus, 
we studied the sensitivity of our results to the location of the 
zones considered. It is therefore important to know how differ- 
ent locations of the R, I, and V zones affect the percentages of 
collisions of planetesimals with the star To test the dependence 
of the results on the location zones, we relaxed the borders of the 
locations by considering 2 cases. In the first case, we moved the 
border between refractory and intermediate from 0.1 AU to 0.3 
AU. In the second case, we moved the border between interme- 
diate and volatile from 1 .56 AU to 1 .2 AU. 



Our results suggest that there are significantly change in the 
percentage values for the CoRot-2 and CoRot-3 systems (see 
second and third columns of Fig. |4|. In these cases, it can be 
seen that a combination of higher and lower values of the ini- 
tial distribution borders R and I (or I and V) results in an abrupt 
transition between refractory and intermediate elements (or in- 
termediate and volatile elements) occurring at the positions that 
we defined as the borders. Thus, a reasonable fraction of parti- 
cles near the borders change from "intermediate" to "refractory" 
elements (or from "intermediate" to "volatile" elements). In the 
first case, the percentage of refractory material that fall onto the 
star is increased by a fraction and the intermediate material that 
fall on the star is decreased by that same fraction. The same idea 
is valid for the second case, a fraction of intermediate material 
near the border is classified as volatile material. Then, the per- 
centages of intermediate and volatile material are increased and 
decreased by the same amount, respectively. 



For CoRot-4 system, there is only a change in the percentage 
values when the border between intermediate and volatile mate- 
rial is moved from 1.56 AU to 1.2 AU. There is no change in the 
percentage of refractory and volatile materials that fall onto the 
star when the border between refractory and intermediate zones 
is moved. This is because the mean-motion resonances play the 
main role in the orbital evolution of the planetesimals and their 
locations are unchanged by the shift of the borders. In this sce- 
nario, there is no change in the distribution of the initial orbital 
semi-major axis of the particles, but only a determination of the 
new location of the R, V, and I zones. 



4.4.2. Planet on eccentric orbit 

Another point for our discussion is provided by the planet's ec- 
centricity ep. We recall that the value of ep for any extrasolar 
planet can vary significantly due to uncertainties in the cuiTent 
observational methods of detection. As we can see in Table |9] 
the v alue of ep for CoRo t -3b and CoRot-2b w ere assumed to be 
zero dAlonso et all 120081: iDeleuil et al.Ll2008l) . Thus, it is impor- 
tant to know whethe r differ ent values of ep yield different evo- 
lutions. IWinter et al.l (|2007') presented a numerical study about 
the influence of ep on capture of planetesimals in mean motion 
resonances during the planet migration. They showed that for 
Mpianet - Mjupiter, whcn cp - 0, the planctesimals are caught 
on 3:2, 5:3, or 2:1 resonances. The 2:1 resonance is responsible 
for the majority of impacts on the star. For ep — 0.1, the 5:2 and 
3:1 are also important, but the 2:1 resonance is strong enough 
to dominate the orbital evolution of the planetesimal. In the case 
ep = 0.2, the 4:1 resonance is important, while the 3:2 and 5:3 
resonances contribute very little to the evolution of the planetesi- 
mals. Other resonances also become important for the cases ep - 
0.3, 0.4, and 0.5. Therefore, the orbital evolution of the disk of 
planetesimals is affected by several resonances but not domi- 
nated by any particular one. In these cases, the planetesimals 
with semi-major axis larger than 0.65 AU are spread without be- 
ing captured in resonance. Until now, we have assumed that the 
CoRot systems have ep - 0. A compar ison of Fig. |4] with the 
results of Fig. 6 from lWinter et all (l2007h for e p = and t = 10^ 
yr show that, in general, the percentage of pl anetesimals that fal l 
on the star is between about of 70 and 80 %. IWinteret al.1 (|2007|) 
demonstrated that this result is a direct consequence of the value 
of the migration speed. 

We also studied the evolution of the system by considering 
the planet to be on an eccentric orbit. In Fig. |6] we present the 
results for ep - 0.1, 0.3, and 0.5, respectively. As expected for a 
circular orbit (see first column of Fig.lUi the percentage of colli- 
sions with the star decreases (from about 10% to 35%) when the 
planet's eccentricity is increased . In general, we note that the 
percentage of intermediate and volatile materials accreted de- 
crease. For the CoRot-3 system, the integrations stopped when 
the last planetesimal was removed from the system. The planet 
did not complete its migration at its current position. The same 
occurred for the planet from CoRot-2 system when ep - 0.5. In 
both systems, the percentage values of refractory material do not 
differ from those for ep - 0. In the CoRot-4 case, some planetes- 
imals had survived when the integrations finished, but the planet 
arrived at its current position, a^f = 0.09 AU. In this case, we 
note an increase in the percentage of refractory material that fall 
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on the star. Planetesimals in the refractory zone are captured in 
mean motion resonances and some of them fall onto the star. 

A planet on an eccentric orbit can significantly increase the 
eccentricities of the planetesimals because of its gravitational 
perturbation force. When planetesimals reach very high eccen- 
tricities they are then removed from the s ystem by collision or 
ejection. As shown in IWinter et al.l (|2007), a fraction of plan- 
etesimals are spread without being captured in resonances. On 
the other hand, other planetesimals are captured in several reso- 
nances (2:1, 3:1, 7:2, 5:1, 6:1, and 5:2). Planetesimals in higher 
resonances are not so well protected against close approaches 
with the planet and a fraction of them are ejected from the sys- 
tem. 



4.4.3. Warm planetesimals disk 

Here we consider a disk of planetesimals with eccentricities ini- 
tially distributed between and 0. 1 . The results are presented in 
Fig. I?] We note that the results are similar to those of a disk of 
planetesimal initially in a circular orbit about a planet with an ec- 
centricity ofO.l (first column of Fig. |6]l. The differences between 
these results and those from the first column of Fig. |4]are due 
to the amount of refractory material that fall on the star. Here, 
only the CoRot-4 system has a percentage value of refractory 
material similar to those for ep =0.0. In the present case, plan- 
etesimals in the refractory zone are ejected from the system. In 
the CoRot-2 and CoRot-3 systems, there is also a decrease in the 
amount of refractory material that falls onto the star, which was 
also ejected from the system. That is because the eccentricity (of 
the planetesimal or the planet) affects the capture in mean motion 
reson ance between the pl anetesimal and the planet. We learned 
from IWinter et alj (l2007h that the 2:1 resonance is responsible 
for the majority of impacts onto the star. According to OuillenI 
(|2006|) . a body is captured into the 2:1 resonance when its ini- 
tial eccentricity is below a certain limit given by en,,, ~ 1.5//'^^. 
This is about 0.126, 0.220, and 0.371 for CoRoT-4, CoRoT-2, 
and CoRoT-3 systems, respectively. 
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Fig. 8. (a)- Relative abundances distribution of elements for the 
CoRoT stars as a function of their condensation temperatures 
(separated), (b)- Superposed distributions: CoRoT-2 (cyan trian- 
gles), CoRoT-4 (blue circles), CoRoT-3 (magenta squares) and 
CoRoT-5 (green stars). 



5. The link between migration and stellar 
abundances 

We now discuss some assumptions and physical conditions that 
are the basis of the link between the consequences of the type 
of planetary migration considered and the distributions of the 
elemental abundances relative to the condensation temperatures 
Tc- To test this link, we present in Fig. [8] the abundances as a 
function of Tc for the CoRoT stars. 

At a certain phase of the disk evolution, we assume that the 
disk is homogeneously filled with planetesimals between a min- 
imal typical distance of ~ 0.03 AU and a maximal one of ~ 4.5 
AU. We also assume that a giant planet is formed at ~ 5 AU, 
before migration. An age of the system (star + disk) of 20 to 
30 M yr is critical for the efficiency of this process dWinter et al.L 
I2007h . At this age, when the disk is already devoid of its gas 
component, a low mass central star will indeed reduce, for the 
first time, its external convective layer to its minimum configu- 
ration (see for instance iFord et al.l (|l99^). In this situation, any 
metallic content of the planetesimals falling onto the star will 
probably remain in the photosphere of the star contributing to its 
metallic enrichment. This may also be true for ages older than 
30 Myr, but not younger than 20 Myr because any enrichment 
would then be drastically reduced due to the larger convective 
layer at these early epochs. 



Clearly, the stellar enrichment depends on the quantity of in- 
fall material. If this mass were of the order of 35 Earth masses of 
planetesimals consisting essentially of iron, this induced migra- 
tion mechanism could furnish an iron enrichment of the order of 
1.8 dex, as observed in stars hosting hot Jupiters (IWinter et all 
:2007.) . But, how would this affect the infall of a collection of 
rocky materials of varying richnesses of refractory, intermedi- 
ate, and volatile elements? First, we must consider that in gen- 
eral, stars apparently without detected planets exhibit flat distri- 
buti ons (zero slope) of the element abundances with respect to 
Tc (lEcuvillon et al.L l200 6l). Then, in principle, any differential 
stellar accretion will produce a positive slope of [X/H] versus Tc 
if the volatile contribution is smaller than the refractory contri- 
bution. In contrast, negative slopes should appear if the volatile 
contribution is far more important. 

Observationally speaking, also in the case of SWF, flat dis- 
tributions are universally found (Ecuvillon et al., 2006). It is 
mainly for this reason that zero gradient distributions are, in 
the literature, synonymou s with the absence o f a self-enrichment 
mechanism. Nevertheless, ISmithetalJ(l200l found, using a col- 
lection of abundance data from the literature, that about six 
stars exhibit p o sitive slopes. Other author s (Gonzalez, 2006a43; 
Sad akane et an. l2002l:lmeda et alillooTI) found no dependence 
of the abundances on Tc- 
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Fig. 7. Same as the first column of Fig. |4l but for planetesimals initially distributed at eccentricities between and 0. 1 . 



Applying our differential accretion mechanism methodology 
to three CoRoT systems (2, 3, and 4) considered in this paper, we 
found surprisingly that in all cases the volatile and intermediate 
contributions were very similar and much larger than the refrac- 
tory ones (Fig. HJi. Here, we present the results fo r a relatively 
slow p lanetary migration of 10^ years. As shown in I Winter et al.l 
(l2007h . for shorter migration times, the accretion mechanism is 



less efficient. We show this in Fig.|5]only for the case of CoRoT- 
4, which involves in the infall of intermediate material only. 

At the top of Fig. |8] we have separated by aiTows along the 
Tc scale, the three V, 1, and R domains. In general, we can in- 
fer that similar V and 1 contributions represent a flat distribution 
of the elements in these two zones. In case of significant stellar 
accretion, the flat distribution will only be shifted upward as a 
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whole. In the case of modest or null accretion, the original and 
flat distribution will be maintained. This is in fact, the case of 
CoRoT systems 2 and 3, where migration, at least in these con- 
ditions is not realized (see Sect. 6). Only in the case of CoRoT- 
4 may differentiate accretion can contribute approximately to a 
mild metal enhancement, but only for rapid migration (Fig.|5]). 

We found that the accretion contribution of highly refractory 
elements such as Ca, for example, is small compared to those 
of I and V elements. In this case, what would be the meaning 
if, in case of a net metallic enhancement, we were to observe 
that the R elements have the same levels of abundances as V 
and I elements? We can indeed infer that R and I elements were 
originally more mixed in the disk as we thought. In other words, 
R elements are contained in I and maybe only the two zones V 
and I are the more realistic representatives ones. 

We can interpret these results in another way. At the end of 
Sect. 4.1, we discussed some effects of the mixing processes 
during the primeval stages of a sol ar/stellar nebula. Th e de- 
tection of silicate bands in comets (iHanner et a n 11994 indi- 
cates that crystalline forms can exist in these outer disk bodies 
(Bouwman et al., 2000; Wooden et al., 2000). Since these crys- 
tals can only be formed at the relatively high disk temperatures 
existing near the star, it has been proposed that a large-scale mix- 
ing occuiTed during the early disk epochs, transporting material 
from the inner to the outer regions of the disk. Mechanisms such 
as turbulence and large-scale meridional flows in regions close 
to the disk midplane can ac t in this way. In con trast, inflow can 
exist above the midplane ( Keller & Gail l2004l) . Whether these 
mechanisms can be changed by the pr esence of magnetohydro - 
dynamical turbulence was studied by ICarballido et alj (l2005h . 
Independently, we have ultimately identified for our case, two 
different situations depending on the extent of these mixing pro- 
cesses. On the one hand, for a partial mixing we maintain the 
R, I, and V zones or only V and I. On the other hand, we ob- 
tain a single large intermediate zone disk when the early disk is 
completely and totally mixed. 

6. Results and conclusions 

In Fig. [8] we present one of our most important result: the rel- 
ative distribution of the abundances of refractory, intermediate, 
and volatile elements as a function of their condensation temper- 
atures for four CoRoT stars. A flat distribution is found for all 
of them, even if CoRoT-3 presents a slight tendency to display 
a positive slope. Clearly, a more reliable determination of this 
slope would be obtained if a complete collection of abundances 
of the highly volatile CNO elements was at hand. In any case, 
we also note the importance to dispose of measurements of the 
volatile S and Zn elements. 

When considering the metallic abundances of CoRoT-2 and 
CoRoT-4, derived in this work, we found that the metallic- 
ity of CoRoT-4 is somewhat higher than that proposed by 
iBouchv et al] (I200 8'). hence confirming a mild metallic excess. 
We compare the stellar parameters with those in the literature 
for CoRoT-2 and CoRoT-4 in Table |2] 

CoRoT-2 star deserves particular attention because of its sig- 
nificant stellar activity and its youth indicators represented, re- 
spectively, by its Ca II H and K lines (Fig. [3]) and its Li reso- 
nance line (Fig. [U, c ontrasting sha r ply w ith those of the older 
1 Gyr star CoRoT-4 dMoutouet aliboOSl) . which does not ex- 
hibit these features. Its age can be determined from the equiva- 
lent width (EW) and abundance of Li by using diagrams of the 
Li values as a function of the T^ff values. By considering its mea- 
sured EW of 124 mA we found, using the diagrams presented 



bv lTorres et all (I2006ll2008h . that this star falls in the lower dis- 
persion limit of the Li abundances found for the members of 
the Pleiades cluster, with an age of 119 + 20 Myr dOrtega et all 
2007). By using its Li abundance foun d here of 2.6, we a rrived 
at the same position in the diagrams of iTorres et al.l (l2006i) and 
Ida Silva et al. I (l2009h . From all these matches, we infer an age 
of 120 Myr for CoRoT-2, from which Bou chv et al.1 (1200 8*) pro- 
posed an age of < 500 Myr. This new age implies that of CoRoT- 
2 is one of the youngest known two stars with planets, the other 
being HD 7 0573, a debris- disk type star of a similar age and Li 
abundance dSetiawan et al.L [2007). The T^g values of these two 
solar-type stars are in a short temperature interval, where old 
solar-type stars with planets exhibit a pec uliar Li depletion wit h 
respect to similar stars without planets dlsraelian et all l2009l) . 
CoRoT-2 and HD 70573 are unaffected by this Li depletion due 
to their youth (see also Pinsonneault (2009). 

We performed simulations of planetary migration into an in- 
ternal disk composed of differentiated particles (planetesimals) 
with refractory, intermediate, and volatile properties, and applied 
these simulations to three CoRoT systems. The accreted material 
contains large and similar contributions of / and V particles and a 
very small contribution of pure refractory elements, as shown in 
Fig. |4] In other words, accretion is mainly "cool" and "warm" 
and not "hot" as largely mentioned in the literature. Because 
zone / exhibits a mixture of chemical properties, mainly contain- 
ing common elements such as Fe, both intermediate and volatile 
contributions produce a flat distribution of [X/H] versus Tq, as 
for the observed three CoRoT systems considered here. 

When there is a complete primeval mixing of elements in the 
disk, a single intermediate zone will be maintained producing a 
flat distribution of abundances as a function of Tq- However, if 
the mixing is incomplete or partial, some R and V rocky mate- 
rial will be maintained near and far from the star, respectively. 
An extended I zone (see top of Fig. [8]) will control a large part of 
the slope. In that case, the agreement of the model with obser- 
vations will depend mainly on the extended distribution of the 
element abundances with Tc < 304 K. Nevertheless, this region, 
containing volatile and highly volatile elements such as CNO 
and noble gases, in general, requires spectra of higher S/N than 
our own to measure the abundances. By considering our O abun- 
dances for CoRoT 2 and 4 and the C abunda nces for CoRoT 3 b y 
iDeleuil et alj d2008h and that for CoRoT 5 dRauer etaUl2009l) . 
we obtained one collection of quite flat gradients as presented 
in Fig. [8] We note that CoRoT 3, for which the gradient is the 
least flat, is a peculiar case because its "planet" is probably a 
brown dwarf star. We conclude that to obtain a reliable collec- 
tion of highly volatile element abundances of CNO and noble 
gases, data of even higher S/N are required. 

In any case, independently of whether accretion is important 
or not, depending on the accreting mass, the main result of this 
work, is that, for the CoRoT systems considered, flat distribu- 
tions that appear to be observationally the rule, do not represent 
the absence of a self-enrichment mechanism as sometimes men- 
tioned in the literature. In contrast, we show here that flat distri- 
butions of the elements abundances as function of Tq could be a 
natural result of accretion. 

Does the discussed type of migration represent true situa- 
tions for the CoRoT systems 2, 3, and 4? This is not the case 
for, at least systems 2 and 3, where unrealistically high disk 
masses are necessary to bring the pl anets to their observed fi - 
nal distances with respect to the star (Adams & Laughlin, 2003!). 
System 4 can be considered an exception^ A different migration 
of type II in a d isk containing gas (Pap aloizou & Terq uem, 200^ 
lArmitageL 12007.) has to be invoked for CoRoT-2. The CoRoT-3 
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system is particularly difficult and could represent a challenge 
to migration theories. If this planetary body, with an exception- 
ally high mass of 21.66 My is a real planet and not a brown 
dwarf, a type II migration study must first avoid the collision 
of the planet with the star due to the large planetary eccentric 
orbits developed, which is the case, at least, for a ~ 10 M/ 
(iRice et al 1 12008"). These authors mention that if the gas disk 
dissipates quickly, the eccentric orbits of these massive plan ets 
could eventually be tidally circularized dFord & Rasiol 120061) . Is 
this the case of CoRoT-3 system? 

CoRoT-4 system with a planet of 0.72 My, also requires a 
high primeval disk mass of ~ 0.2 - 0.3 M©. Even if it were 
an extreme case, it could be acceptable for a star of 1.2 Mq. 
Nevertheless, our simulations indicate that in this case only a 
very rapid migration of 1000 yr can produce a metal enhance- 
ment as observed of [I/H] ~ 0.1., if occurred in the first 20-30 
Myr, when the stellar conv ective layer for the first time attains 
its minimum configuration (iFord et alill999,) . This could also be 
possible in principle, ev en for the smaller p lanet of mass 0.467 
My related to CoRoT-5 (iRauer et all 12009 1). We also show that 
these results are highly dependent on the model adopted for the 
disk distribution regions in terms of refractory, intermediate, and 
volatile elements and the other parameters considered. 
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